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Abstract —A continuous-time Wiener phase noise channel with 
an integrate-and-dump multi-sample receiver is studied. A lower 
bound to the capacity with an average input power constraint 
is derived, and a high signal-to-noise ratio (SNR) analysis is 
performed. The capacity pre-log depends on the oversampling 
factor, and amplitude and phase modulation do not equally 
contribute to capacity at high SNR. 

I. Introduction 

Instabilities of the oscillators used for up- and down- 
conversion of signals in communication systems give rise to 
the phenomenon known as phase noise. The impairment on 
the system performance can be severe even for high-quality 
oscillators, if the continuous-time waveform is processed by 
long filters at the receiver side. This is the case, for example, 
when the symbol time is very long, as happens when using 
orthogonal frequency division multiplexing. 

Typically, the phase noise generated by oscillators is a 
random process with memory, and this makes the analysis of 
the capacity challenging. The phase noise is usually modeled 
as a Wiener process, as it turns out to be accurate in describing 
the phase noise statistic of certain lasers used in fiber-optic 
communications 12. As the sampled output of the filter 
matched to the transmit filter does not always represent a 
sufficient statistic (2|, 0, oversampling does help in achieving 
higher rates over the continuous-time channel 0-0. 

To simplify the analysis, some works assume a modified 
channel model where the filtered phase noise does not consider 
amplitude fading, and thus derive numerical and analytical 
bounds l7ll— ifTOl . 

The aim of this paper is to give a capacity lower bound 
without any simplifying assumption on the statistic of filtered 
phase noise. Specifically, we extend the existing results for 
amplitude modulation, partly published in Q, and present new 
results for phase modulation. 

Notation: Capital letters denote random variables or random 
processes. The notation X ^ = (X m , X m+ i ,..., X n ) with 
n > m is used for random vectors. With Af(0, a 2 ) we denote 
the probability distribution of a real Gaussian random variable 
with zero mean and variance a 2 . The symbol = means equality 
in distribution. 

Given a complex random variable X, we use the notation 
A’ and ZX to denote the amplitude and the phase of 


X, respectively. The binary operator ® denotes summation 
modulo [—7r,7r). 

The operators E [•], h (•), and /(•;•) denote expectation, 
differential entropy, and mutual information, respectively. 

II. System model 

The output of a continuous-time phase noise channel can 
be written as 

Y(t) = X(t)e j6W + W(t), 0 < t < T (1) 

where j = \J 1, A (■) is the data bearing input waveform, and 
W is a circularly symmetric complex white Gaussian noise. 
The phase process is given by 

0(f) = 0(0) + 7 y/TB{t/T), 0 <t<T, ( 2 ) 

where B(-) is a standard Wiener process, i.e., a process 
characterized by the following properties: 

. B(0) = 0, 

• for any 1 > t > s > 0, B(t) — B(s) ~ Af(0, t — s) 
is independent of the sigma algebra generated by 
{B[u) : u < s}, 

• /i(-) has continuous sample paths. 

One can think of the Wiener phase process as an accumulation 
of white noise: 

0(f) = 0(0) +7 [ B'(r)dT, 0 < f < T, (3) 
Jo 

where B'(-) is a standard white Gaussian noise process. 

A. Signals and Signal Space 

Suppose X(-) is in the set £ 2 [0,T] of finite-energy signals 
in the interval [0, T\. Let {</> m (f )}“ =1 be an orthonormal basis 
of £ 2 [0,T]. We may write 

oo oo 

X(t) Mt), W(t) = W m (4) 

m =1 m= 1 

where 

X m = f x{t)<t> m (t)*at, (5) 

Jo 

x * is the complex conjugate of x, and the {W m }” =1 are 
independent and identically distributed (iid), complex-valued. 


circularly symmetric, Gaussian random variables with zero 
mean and unit variance. 

The projection of the received signal onto the n—th basis 
function is 

Yn= [ Y(t) (j) n (ty dt ( 6 ) 

Jo 

oo 

= £ x 
m— 1 
oo 

= El X m &mn + W n . (8) 

m—1 

The set of equations given by © for n = 1,2,... can be 
interpreted as the output of an infinite-dimensional multiple- 
input multiple-output channel, whose fading channel matrix is 


f e^dt + W n (7) 

'n 


B. Receivers with Finite Time Resolution 


Consider a receiver whose time resolution is limited to A 
seconds, in the sense that every projection must include at 
least a A-second interval. More precisely, we set MLA = T, 
where M is the number of independent symbols transmitted 
in [0,T] and L is the oversampling factor, i.e., the number 
of samples per symbol. The integrate-and-dump receiver with 
resolution time A uses the basis functions 


for m = 1 With the choice ©, the fading 

channel matrix ( 1» is diagonal and the channel’s output for 
n = 1,..., ML is 

7.A 


Y n = X n 1 


e jQ(t) dt + Wn 


(n— 1) A 


= I n e J ' e(( "- 1)A) ^ [ e j(e(t) " e(( ”" 1)A)) df + W, 

X J (n— 1) A 

= X n e jQn ]- f e jn/VKBn( ' t / A) dt + W n (10] 

A Jo 

= X n e j@n [ e J7v / AS„(t) dt + Wn 

Jo 


= X n e j& "F n + W n 


( 11 ) 


where we have used the notation 0„ = 0((n — 1)A) and 
F n = fy e J 7 v / AB„(t) j n m we have usec j ([ 0 ]) property 

B(t/T) - B((n - 1)A/T) = B(t/T - (n - 1)A/T), the 

substitution 

t <— t — (n — 1)/A 
B n {t/T) «- B(t/T - (n - 1)A/T), 


( 12 ) 


and the property s/TB n {t/T ) = \fAB n (t/ A). Finally, in 
step (a) we have used the substitution t <—t/ A. 

Since the oversampling factor is L, we have X k L+i = 
XkL +2 = ... = XkL+L for k = 0,..., M — 1, and we can 
write the model as 


Yn = X ln/L]L e^F n + W n 


for n = 1,,ML. 

The vectors X{ vr/j , F j M L , and Wf* L are independent of 
each other. The variables {X k i J }^L 1 are chosen as iid with 
zero mean and variance E [|X„| 2 ], and the average power 
constraint is 


E 


1 

T 



\X(t)\ 2 dt 


1 


ML 


MLA 

E[\X„ 


E E D Y « 


n— 1 
121 


< v. 


(14) 


Since we set the power spectral density of W to 1, the power 
V is also the SNR, i.e., SNR = V. 

Using ©, the variables (-) j 1L follow a discrete-time Wiener 
process: 


©n — ©n-i + N n -i, n—l,...,ML, (15) 

where the N n ’s are iid Gaussian variables with zero mean and 
variance 7 2 A. The fading variables F n ’s are complex-valued 
and iid, and F n is independent of 0". In other words, F n 
is correlated only to N n , and is independent of the vector 

(ivr 1 ,^). 

Note that for any finite A, or equivalently for any finite 
oversampling factor L, the vector Yj Mi does not represent 
a sufficient statistic for the detection of X given Y in the 
model ©. 


III. Lower bound on capacity 

We compute a lower bound to the capacity of the 
continuous-time Wiener phase noise channel (IT3li- ([T5l i. For 
notational convenience, we use the following indexing for 
i = 1,..., L and k = 1,..., M: 

Y(k-i)L+i = X k 3i+i F'( fc _ 1 ) Il+i + W(k-i)L+i, (16) 

and we group the output samples associated with X k in the 

vector Yi. — 

vector 1 k — 1 (k-i)L+i ■ 

The capacity is defined as 

C ( S N R ) = L sup / (Xf ; Y, M ) (17) 

where the supremum is taken among the distributions of X^ 1 
such that the average power constraint (IT4l) is satisfied. 

The mutual information rate can be lower-bounded as 
follows: 

it I (*1 M ; ■Yf) = if E 7 (**; Y 1 M IE' 1 ) 

/c=l 

1 M 

- ME 7 ( Y "I x i~ l ) + 1 (Yf I x*~\\Xk\) 

1 k=1 

( b ) i yi 

> jj E 7 (|Yfc| 2 ;H Y ^II 2 ) + 7 (^Y fe ; Ytt I X k ~i, \x k \) 

k= 1 

= I (lYtl 2 ; IIYiH 2 ) +/(ZX i; Y^ | X 0 , |*i|) (18) 

s— _ ✓ s— ✓ 


(13) 










where step (a) follows by polar decomposition of X k , step (b) 
holds by a data processing inequality, by reversibility of the 
map x i—>■ x 2 for non-negative reals, and because 1 is 
independent of (X k ,Y k ). Finally, the last equality follows by 
stationarity of the processes. 

A. Amplitude Modulation 

By choosing a specific input distribution that satisfies the 
average power constraint we always get a lower bound on the 
mutual information, so we choose the input distribution as 


By putting together (l23l > and < 125b we obtain 


> ———-F - ln(L 2 /i 2 A 2 + Xu) 

- - (SNR-Var [G] +2 +A 4 ) 


2 ln 


( 26 ) 


In the limit of large time resolution we have 


lim 

A-ro 


Var [G] 
A 3 



(27) 


P\x k \ 2 ( x ) 


I exp (— ) x > A * 

0 elsewhere 


(19) 


where A = SNR A — A 4 > 0 with t > 0. Note that with this 
choice the average power constraint is satisfied with equality, 
i.e., E [\X k \ 2 ] = SNRA. 

Similar to the method used in a, we give here a lower 
bound to the first term on the right hand side (RHS) of ( ITSl i 
in the form 


I\\ > E[-lnq v (V)]-E[-lnq vllXll2 (V\ l^l 2 )] (20) 

where V = 11Y 1 11 2 and 


Now we let the time resolution grow as a power of the SNR, 
i.e.. A -1 = [SNR“], and the parameter u = pA~@, with 
p > 0. By using (E71 > into (E6l) . in order to find a tight bound 
in the interval 1/3 < a < 1 we need to satisfy the condition^ 
a < l/(f +1) and /3 > 1. The tightest bound is obtained with 
/3 = 1 and p = 4: 

sn'Soo^" - ^ ln ( SNR )} ^ ^^ ln ( 47re )' ( 28 ) 

For 0 < a < 1/3 we need to satisfy the conditions a < 
l/(i+l) and a > l/(/3+2), and the tightest bound is obtained 
by choosing /? = a -1 — 2 and p = 27 2 / 45 : 


qv(v) 



P\x 1 p{x)q V \ \ Xl \*{v\x) dx. 


( 21 ) 


lim 

SNR— Foo 



^ hr (SNR)} 



/ 27T7 2 e\ 

J ■ 


(29) 


Specifically, we choose the auxiliary channel distribution as 


1 


Qv \\ Xl \*{v\x) = ^^exp 


(x-L(l + xE[G])) g 


where G = ||Fi|| 2 /L and v > 0, for which we havf] 

ln<3V| |.Y 1 | 2 (1 / '| 

L ( E [|Xr| 2 ] 


( 22 ) 


E [-lnq vllXll i(V\ lAil 2 )] = ± ln (nu) + l E [lnd^l 2 )] 


-Var [G] + 2E [G] + E 


1 


Ai| 


< i ln (ttuX) + — r E-—|- — (SNR • Var [G] + 2 + A 4 ) (23) 
L zA v 

where the inequality is due to E [G] < 1, E [|Xi| 2 ] < SNRA, 
the bound E [|Xfc|~ 2 ] < A 4 which follows from the support 
of |Xfc| 2 , and 


r°° i 

E [In | Ai| 2 ] = / -exp 
J A-* A 


- A“ 


ln A - 


< ln A + 


JA-t/x 

A" 4 


exp — it — 


ln(x) dx 

A" 4 


A 


A ' 


ln(u) dit 


(24) 


By substituting (l22l > and ( 1 1 9b into (EH . and by following 
similar steps to those of 0, we get 

E[-ln< 2 y(V)] > i ln(L 2 p 2 X 2 + Xu). (25) 


B. Phase Modulation 

The second term in the RHS of (ITSl) can be lower-bounded 
as follows 


(a) 


I (^Xj_; Yq | A 0 , |Ai|) > I(ZX 1 ;$\X 0 ,\X 1 \) 

> E [- In|( < T > |JAo, |Xi|)] - E - lng^^^lAg) 

(30) 


where step (a) is due to a data processing inequality with 

$ = Z(Fi(F 0 e“ J ' ZXo )*) 

= ZXj. © /((Ail F, + W,) 0 Z(\X 0 \FZe jN ° + W 0 *), 

(31) 


and the last inequality follows by choosing the auxiliary 
channel 




■|*o) = 


exp(£ cos(£> — Zxi)) 

2Wo(0 


(32) 


where Jo(-) is the zero-th order modified Bessel function of 
the first kind, and £ is a positive real number. Since we assume 
an uniform input phase distribution, the output distribution is 
also uniform: 


fe|A 0 ,|Y 1 |(^|a;o,|a:i|) 



1 

27r 

(33) 


'Details are provided in the extended version of the paper. 


2 Details are provided in the extended version of the paper. 





















Using (l32l >. the second term in the RHS of ( 13 01 1 can be upper- 
bounded as follows for any A < A < cxd: 


=ln(27r/ 0 (C)) - [cos($ - ZX x )\ 

< ln(7rV7r) + ^ ln + C P 


= 2 ln { 2 ^ep) 


(34) 


where the inequality is due to Io(C) < yA "/ 2 ' e< */V( derived 
in flTl Lemma 2], and from the result of Appendix lAl with 

P = 1-E [F 0 e~i N °] E [F!] + 2e- 3 A A / 8 E [IXil” 2 ] K k (35) 

where K A > 1 is a finite numbefl The last step in m is 
obtained by choosing ( = ( 2 p) _1 . 

In the limit of large time resolution we have 


lim 

A—>o 


i- 

lA 


- 2K a A 




(36) 


where the inequality follows from the bound E [| X± | 2 < A*. 
Choosing t = 1 and putting together (l30l > and < 133 b -( 136b we 
get 


Jim {Iz + ^ln(A)} > 7 In 


A—X) 


1 




7te(7 2 + 3iC A ) 


(37) 


and letting the time resolution grow as a power of the SNR, 
i.e., A -1 = [SNR“], for 0 < a < 1/2 we have 


lim 

SNR—^oo 


{/z - % In 


(SNR)} 


^ ln 


7re(7 2 + 3 K A ) 


(38) 


IV. Discussion 


As a byproduct of <f28l) . ( 1291 . and (l38l >. a lower bound to 
the capacity pre-log is 


C (SNR) 

lim —- 

SNR—s-oo ln(SNR) 


> 


2a 0 < a < 1/3 

(1 + a)/2 1/3 < a < 1/2 (39) 

3/4 1/2 < a < 1. 


Figure |T| shows the lower bounds on the capacity pre-log 
versus the parameter a, as reported by ( l39l >. The contributions 
of amplitude and phase modulation are also shown separately: 
Amplitude modulation reaches full degrees of freedom by 
sampling more than vSNR samples per symbol, while phase 
modulation achieves at least half of the available degrees of 
freedom by using a time resolution that scales as 1/v/SNR. 

The input distribution that achieves the capacity lower 
bound is uniform in phase and the square amplitude is dis¬ 
tributed as a shifted exponential ( IT9l >. The statistic used for 
detecting \Xk\ is ||Yfc||, and the one used for detecting ZX^ 
is Z {Y[k- t)Le _j 7 ' Yfc -i) }■ 



Fig. 1 . Capacity pre-log lower bounds as a function of a at high SNR. The 
oversampling factor L is L = [SNR“], 


V. Conclusions 

We have derived a lower bound to the capacity of 
continuous-time Wiener phase noise channels with an average 
transmit power constraint. As a byproduct, we have obtained a 
lower bound to the capacity pre-log at high SNR that depends 
on the growth rate of the oversampling factor used at the 
receiver. If the oversampling factor grows proportionally to 
SNR a , then a capacity pre-log as high as that reported in ( l39l > 
can be achieved. 


Appendix A 

A LOWER BOUND TO E [cOs($ - ZX i)] 
The expectation can be simplified as follows: 


E [cos(4> — ZX i)] 

( => E [cos(Z(|X 1 |J ’ 1 + Wi) - Z(\X 0 \F 0 e-3 N ° + Wo))] 

( = } E [cos(Z(|X 1 |i 7 ’ 1 + Wi))] E [cos(Z(|A 0 |F 0 e-^ + W 0 ))] 
+ E [sintZflXylFi + Wi))] E [sin(Z(|Xo|F 0 e- jAro + W 0 ))] 

= E [costZflXilF! + Wi))] E [cos(Z(|X 0 |Foe _jJV ° + W 0 ))] 

(40) 


where step (a) is due to <f3Tb . step ( 6 ) to the addition formula 
for cosine and independence of random variables, and the last 
step follows because E [sin(Z(|Xi|Fi + W\))\ = 0 as F\ and 
W\ have symmetric pdfs with respect to the real axis. 

The first expectation on the RHS of (l40t can be written as 


5R{e- 


jA{\X 1 \F 1 +W 1 ) 


sft{ e i ZFl e jZ ( lXlFl \ +Wl) } 


= E [cos(ZF 1 )cos(Z(|X 1 F 1 | + W x ))] 

(41) 


3 For example, choosing 7 2 A = 0.01 gives K A = 8 . 1353 . See the where the first ste P is due t0 the circular symmetry of Wi, 
extended version of the paper for a detailed derivation. and the second Step because of the symmetric pdfs of F\ and 























W\. A lower bound to SB is given by 


( b ) 

> E 


^ e jAQX 1 \F 1 +W 1 )^ ( > } E [sr{f 1 }cos(Z(|XiFi| + Wi))] 

1 


(c) 

> E 


(d) 

> E 


> 


5i{Fi}l0ft{Fi} > 0) ^1 - 
+ E [5ft{Fi}l(K{Fi} < 0)] 
3?{Fi} - l(5R{Fi} > 0) 

1 


l-MI 


|*i F, | 2 


K{Fi}- 




7 2 A 


(1 - e ~ l2A/2 ] - E 

1 

V ) 

U*iN 


A', 


(42) 


where step (a) holds because |Fi| < 1, (6) follows by 
cos(at) < 1 and b>0 


E[cos(Z(p + Wi))] > 1- -4, p > 0, (43) 


step (c) because 5i{Fi} < 1, step (d) is obtained by subtract¬ 
ing E [l(5ft{Fi} < 0) |AiFi| -2 ], and the final inequality uses 
E [| Ai| -2 ] < K A for a finite suitable A0. 

Following an analogous derivation used for finding (l42l> . for 
the second factor on the RHS of (l40t we have 


»{ 


o jZ(\X 0 \F 0 e- jN °+W 0 ) 


> E 


3£{F a 


,-jNo 




1 


> 


2-7T 

7 2 A 


erf 


7 2 A 


-37 2 A/8 _ E 


U*ol 2 J 


|*oF 0 | 2 
K a (44) 


where erf(-) is the error function, and the closed form for 
E [Foe _jAr °] is provided in Appendix iBl Using (l42l) and < |44]> 
into ffl, with 5R{E [Fi]} < 1, 5J{E [F 0 e~ jNo ]} < e~ 3 ^ A / 8 , 
E [|Xo| -2 ] > 0, and E [|Fq|~ 2 ] > 0, the final result is 


E [cos($ - ZXi)] > E [F 0 e~ jN °] E [Fi] 


- 2e- 3 ^ A / 8 E 


l*i I 2 


K 


A- 


(45) 


Appendix B 

Evaluation of E [ F 0 e~ jN °] 

Knowing that Nq = a B(t) dr with er = 7 \/A, we 
compute 


Var [crB(t) — Nq] 


<j 2 \lar [B(t)] -f Var [Nq] — 2aE [B(t)No] 
a 2 (t + 1) - 2ct 2 [ E [B(t)B(r)] dr 

J o 

a 2 {t 2 -t+ 1) (46) 


where the last step follows from the property of Wiener 
processes E [B{t)B{r)\ = min{f,r}. Thus we have 


E [F 0 e- jJVo ] 



e j(rrB(t)-N 0 ) 


dt 


Var[rTS(i)-JV 0 ]/2 fa 



(47) 


where in step (a) we used the characteristic function of a 
Gaussian random variable, and in the last step we used (l46t . 
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Appendix C 
Auxiliary channel 

Choose the auxiliary channel distribution 

, . . 1 ( (v - L(1 + x/z)) 2 \ 

= - ~ x -1 (48) 


where 


and 


G = 


|Fi| 


L ’ 

Squaring the output statistic gives 


j=i 


C 2 = |Xi| 4 L 2 G 2 + 4\X x \ d LGZ x + 2\X 1 \ 2 (2 Zf + LGZ 0 ) 


+ 4\X x \ZqZ x + Zq, 


thus 


H^i| 2 

T 2 AT 

= — E [|X!| 2 ] E [(G - y) 2 ] + —E [|X!|] E [Zi(G - M )] 


+ -(2E[Z 2 ] +LE[(Z 0 -L)(G-/x)]) 

4 r i i i 

E[Z!(Z 0 -L)] + -E 


H—E 

V 




i 


l*i| 2 


The conditional entropy is 

E [- lngy| i.XipCFI |Xl| 2 )] 

= \ H™) + |e [ln(|Xi| 2 )] + E 

= iln(TO/) + ^E [ln(|ATi| 2 )] 


(c-ca + i^iv)) 2 

h^i 2 


i /e[|Xi| 2 ] 


Var [G] + 2E [G] 


U*il 2 J 


V = IIYiH 2 = \X x \ 2 LG + Z 0 + 2\X x \Z x (49) 


1 A _t 


+ 


SNR - Var[G] +2E[G] 


l*i| 2 J 


(54) 


Z 0 = ||W 1 || 2 , Z!=J2^ w *}- 

(50) 


where in the last inequality we have used the bound (l24l >. 
The output distribution is 

pOO 

Qv(v) = p\ Xl \ 2 (x)qv\\ Xl \ 2 {v\x)dx 
Jo 


1 


/ 0 

f°° 1 / x-A~^ 

/ t ex P (-r~ , , 

/ A -t A \ A y yjKVX 

exp (JjLzW±sAl ]dx 


(51) 


< exp 


1 


(V)i r xp (-!)vsi 


(F-ca + i^iv)) 2 

= -C 2 - 2CL(1 + lATil 2 ^) + L 2 (l + IAGIV) 2 
= |Ai| 4 L 2 G 2 + 4|Ai| 3 LGZi + 2|Xl| 2 (2Z 2 + LGZ 0 ) 

+ 4\Xi\ZqZi + Zq + L 2 (l + \X x \~ y) 2 
- 2(\X 1 \ 2 LG + Z 0 + 2\X x \Z x )L(l + \X x \ 2 y) 

= |Xi| 4 L 2 (G - y) 2 + 2|A 1 | 2 (2Z 2 + L(Z 0 - L)(G - y)) 

+ 4\X x \ 3 LZ x (G - y) + 4\X x \Z x (Z 0 - L) + (Z 0 - L) 2 . 

(52) 

Taking expectations gives 


(V-L(l + \X x \ 2 y)f 


exp ( - 

'A~* 


X J Jo A 

(V - L( 1 + xy)) 2 


dx 


= exp 

exp 

= exp 

exp 


1 


A J J q LyX 
{v- L-y)) 2 
vy/{Ly) 


exp 


LyX) ^/niAy/(Ly) 


d y 


A- 


1 


A / \JLyX(LyX + z//(L/z)) 
/ 2L/z 


\ ^ 


1 — T — \v — L L /1 + 


L 2 y 2 X 


(55) 


and the entropy of the output of the auxiliary channel is 


A 


1 


E [- lngv(C)] =-— + - ln(LyX(LyX + v/{Ly))) 


2 Ly 


E[V-L\-E[\V-L\\Jl + 


L 2 y 2 X 


(a) A _i 
> 


A ■- 2 ln(LpA(LpA + v/{Ly))) 


+ —E [V-L] 


E [(2 0 - T) 2 ] 

1 


> — 


v 

A-‘ 


1 + 


L 2 y 2 A 


- 1 


+ -ln(L 2 p 2 A 2 + Ai/). 
A z 


(56) 


= VE[^| 2 ]E[(G-,) 2 ] + vE[G]^E L | Xi | 2 j 


( 53 ) 


where step (a) holds because E [| ■ |] > E [•], and the last 
inequality holds because of 

E [V-L\ = E[\X x \ 2 LG + Z 0 +2 \X x \Z x -L] 

= E [\X x \ 2 \\F x \\ 2 ] >0. 


which is minimized by choosing y = E [G]. 


(57) 































































The lower bound to the mutual information rate for the In order to have a tight bound, we need to satisfy the 
amplitude modulation is 


constraints 


I\\ > E[-ln q v (V)} - E [— In q v \ | Xl | 2 (V| |Xt| 2 )] 

ra-* I 1 

>-:-1— ln(L 2 /i 2 A 2 + Av) -In (nisA) 


2A 

-- (sNR-Var[G] + 2E[G] + E 


Ll^i 


• (58) 


1 — a(t + 1 ) > 0 
a(l-/3) + l-3a<0 
a(l -0) < 0 
a(l — 0) — at < 0 
1 + a(l — /9) > a{t — 0 + 2) 


(62) 


Using A = SNRA — A 4 and choosing v = pA and A = that reduce to 
L- 1 = SNR" Q gives 


In > - 


1 


1 


In (SNR 2a (E[G]) 2 


2 SNR 1- “ (t+1) — 1 2 

(SNR 1- “ - SNR Qt ) 2 +p(SNR 1 " a(1 “ /3) - SNR“ (t+/3) )) 


a < l/(f + 1) 
a(l-/3) + l-3a<0 
a{l-0) < 0 . 


(63) 




- - In \npSNR 1+a V~ 1 '> 


1- 


1 


snr 1_ “( 4+1) 


- p~ l SNR a(1_/3) (SNR ■ Var [G] + 2E [G] + SNR' at ) . 


Next we consider the two cases 1/3 < a < 1 and 0 < a < 
1/3. 

If a > 1/3, i.e., 1 — 3a < 0, then we have to satisfy 


For large SNR we have 


lim E [G] = 1, limAH = 2! 

A— >0 1 J A— >0 A 3 45 


(59) 


(60) 


a < l/(f + 1) 

0 > 1 . 


(64) 


to get 


that gives 

„ f 3 1 

lim /(|Ai| 2 ;U)> lim < — ---——- 

SNR—^oo V ' SNR—>oo ( 2 SN R 1_ “( t+1 ) — 1 

+ i In ((SNR - SNR“ (t+1) ) 2 + p{ 'SNR 1_a(1-/s > - SNR^+^M 


SN 1 i°A{ J (' X -' a ^)- 1 + Q 2~ ii> i “( SNR >} 


1 , . , 2 
> ——In [np) - 
2 P 


(65) 


- 4 l n ( 7 r j oSNR 1+a( ^ _1) 


1 


SNR 1 ~“ (t+1) 


nd the tightest bound is obtained by choosing 0 = 1 and 

P = 4: 


- p _1 SNR“^ 1- ^ ( SNR 1_3a • — + 2 + SNR- 
1 45 


{-T 0^1= - 1 to (5NR)} > -i ‘ n ( 4 ^e). 


= lim i - - 


1 


SNR—^oo 2 SNR 1_Q(t+1) - 1 

„ , ,, , m\ If a < 1/3, i.e., 1 — 3a > 0, then we have to satisfy 

? a(t+i)^2 SNR 1- — <;mp n ( 4 +p) \ ' 1 


1, / (SNR-SNR^ 14 ) 2 


SNR 1+ “^- 1) 


2 ln 


- g ln I 71 P ( 1 - 


p- 


— SNR° 


SNR 1_ “ (t+1) 

( 2 
SNR 1_3 “ ■ — + 2 + 
45 


SNR 


SNR- 


l+a(/3-l) 


a < l/(t+l) 
a > 1 /(0 + 2). 


( 66 ) 


(67) 


= lim 


SNR—^oo 2SNR 1 ~ Q(t+1) - 1 

+ i ln (SNR 1+a(1_/3) - 2SNR“ ( ‘“ /3+2) + SNR a(2t “^ +3) " 1 
P 


By choosing 0 = a 1 — 2 and p = 27 2 /45 we get 
snr-Ioo { 4 A. I 2 ;V) - f 1lo ((5NR)} > a :i» ) . 


( 68 ) 


+P~ 


1 


SNR 


l-a(t+l) 


— — In ( 7T/3 (1 — 


SNR 


1—a({+l) 


- 1 QMR“( 1 -/ 3 ) 


-p- L sm 


SNR 1_3a • —+2 + SNR"“* 
45 


( 61 ) 


Appendix D 

A LOWER BOUND TO E [cOs(Z(p + W))} 
The pdf of 4/ = Z(p + W) is fl2l 


P^Up) = —e p + _ pcosUp)e p sin ^ erfc(—pcos(^O) 

27r V47T 

( 69 ) 























where erfc : a; i—>- 1 — erf(x) is the complementary error 
function. A lower bound to E [cos(’I')] if]: 


E [cos(H/)] == 2 I cos(tp)p^('i/j) d’lp 

( b ) r p 


f‘TX 

/ —= cos 2 ('i/j)e~ p sin ^ erfc (—p cos(-0)) d0 

J 0 V 77 

(c) /*7r/2 ^ 

^ / -t—cos 2 ('ip)e- p sin ^ erfc(—pcos(^)) 

Jo V 77 


> 


(d) W2 „ 
> 1 1 


[ 4 = cos 2 (ip)e~ p2 sin 2 W (2 - e ~p 2 ™ 2 W) 

Jo V 77 


t /2 


cos 2 {iP)e~ p2 sin2w cty - ^-pe~ p2 


(70) 


where (a) follows by symmetry, (&) follows by using (l69l > 
and f Q cos( , 0 ) dt/’ = 0 , (c) holds because the integrand 
is non-negative over the interval [7r 2 ,7r], ( d ) holds because 

erfc(/ocos(V0) > 2 -e~ p2cos2 W0 and cos 2 (t/0e _p2 sin2W > 0 , 

and finally the last step follows by direct integration. 

We bound the integral in (TTOt as follows 



^=cos 2 {^)e~ p2 sin2 W di/> 
V 77 






1 

V 


(71) 


where inequality (a) follows from cos 2 (t/>)e _p Sln W > (1 — 
t/j 2 )e~ p ^ , and the last step holds because erf(-) < 1 and the 
last term is non-negative. Substituting back into (f70t yields 


E [cos(*)] > erf (y) 


1 

v 




-pe 


(a) 




(&) 

> 1 - 




_4_1_ _ f_3_\ 

7 r 2 p 2 e 2p 2 4p 2 \2e/ 

1 


3/2 


(72) 


where step (a) is due to erfc(x) < e 3:2 , and inequality ( b ) fol¬ 
lows from p 3 e~ p < (3/(2e )) 3 / 2 and e -7r p / 4 < 4 /(- 7 r 2 p 2 e). 


Appendix E 

An upper bound to E [|.Fi| -2 ] 

Denoting Z = F\ \, we compute an upper bound as follows 


E [Z- 2 ] 


= lim 
<40 



-^Pz{x) dx 
x z 


(a) 


lim 

<54.0 


Pr(Z < 6) 
P 



Pr(Z < x) dxl 
x 6 J 


(b) ,00 9 

< / — E [1(Z < x)] dx 

Jo x 

= [ — E [ 1 (Z < x)] dx + 
Jo x 

-J 0 + h 



2 


E [1(Z < x)] dx 


(73) 


where e is a suitably chosen positive number, step (a) follows 
by integrating by parts, inequality ( b ) holds because the cumu¬ 
lative function is always positive, and the last inequality holds 
by choosing a function g(Z) > 1(Z < x) for the first integral 
(i.e., for 0 < x < e) and the inequality E [1(Z < x)] < 1 for 
the second integral, that can be computed in closed form. 

As for the function g(Z) we choose 


g(Z) = a(l - Z 2 ) (1 - Z 2 Pl ) (1 - Z 2 P2 ) 

= a [l — Z 2 (1 + pi + P 2 ) + Z 4 (pi + P 2 + P 1 P 2 ) — Z 6 pip 2 ] , 

(74) 


a polynomial whose positive roots are { 1 , p x l ^ 2 1 p 2 1,/2 }> and 
with 3 ( 0 ) = a > 1. To guarantee the positivity of g(Z) for 
0 < Z < 1, we can design the roots p 1 and p 2 to be 
greater than 1 , hence p\ and pi less than 1 . 

If the polynomial g(Z) satisfies the condition E [ g(Z)\ = 0, 
then we have a finite bound in (f73l) : 


E M 

Imposing the condition E \g{Z)\ = 0 in (l74l > gives 

— 1 + E [Z 2 ] (1 + Pi) — E [Z 4 ] Pl 
P2 “ —E [Z 2 ] + E [Z 4 ] (1 + pi) - E [Z 6 ] p! • 


(75) 


(76) 


We want P 2 to be positive, for this we distinguish two cases. In 
the first case we have both numerator and denominator of (f76l > 
positive, and this is satisfied if 


E [l — Z 2 ] E [Z 2 (l — Z 2 )] 
E [Z 2 (l — Z 2 )] ’ E[Z 4 (1-Z 2 )] 


>1, (77) 


so this situation is not wanted. The other case is where both 
numerator and denominator are negative, i.e., for 


E [l — Z 2 ] E [Z 2 (l — Z 2 )] | 

E [Z 2 (l — Z 2 )] ’ E [Z 4 (l - Z 2 )] J ' 


(78) 


Moreover, we want P 2 < 1, and imposing this condition 
on (f76t when numerator and denominator are negative means 

E [(1 ~~ Z 2 ) 2 J 
Pl ~ E [Z 2 (l — Z 2 ) 2 ] 


6 The proof is the one proposed in the Ph.D. thesis of H. Ghozlan. 


( 79 ) 











Conditions (f78l) and (l79t can be numerically checked by 
considering that 0 

E [Z 2 ] = ^2 (—1 + e ~ at + at) (80) 


E [Z 4 ] = 


392 


87 

T 


40 

- y ae- 


18 


e~ 4a - 30a + 8 a 2 


(81) 


E [Z 6 ] = 


— 100a 3 e a + 144a 3 — — ae 4a 
25 


11991 


-ae 


1499a-ae 

200 


—9a: 


2123 2„-a , 4 2„-4a 7nn 2 
—-—a e 4—-a e — 792a 

3 15 


(82) 


where a = 7 2 A/2. For example, for a = 1.3 and jy/A = 0.1 
the roots are p 1 1 l,p 2 « 1.0001 > 1 . 

The last thing to check is that 


g(Z) > 1, for 0 < Z < x and all 0 < x < e. (83) 

Studying the convexity of g(Z), we find that the function is 
concave in 0 < Z < e n with 

2 _ Pi + Pi + P\P2 / / Pl+J>2+_PlP2_\ 2 _ 1 + Pl + P2 

5pip 2 V V 5 Pl/°2 / 15P1P2 

(84) 

Moreover, we have g(e i) = 1 with e\ given by Cardano’s 
formula. Condition (l83l > is verified if e < min{ei,e n }, 
because for all 0 < Z < e we can guarantee g(Z) > 1 
thanks to g(0) > 1, g(e) > 1, and concavity for g{Z). For 
example, for a = 1.3 and 7 \/A = 0.1 we have e\ « 0.3506 
and e n ~ 0.5774, so we choose e = 0.3506. Using (1751 1 this 
gives a bound E [|.F| -2 ] < 8.1353 for all jVA < 0.1. 
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